a=load('f16_0_5sphere_0_4offset_mlf.rcs');
b=load('fixed_f16_0_5sphere_0_4offset_PO-MoM_300MHz.dat');
c=load('fixed_f16_0_5sphere_0_4offset_EPA-PO-MoM_300MHz.dat');
tick=[-180:1:180];

figure(1);
plot(tick,a(:,5),'b-');
hold on;
plot(tick,c,'r+-','MarkerSize',3);
xlabel('\theta_s');
ylabel('RCS(dBsm)');
legend('MLFMA','PO-MM with EPA');
title('VV polarization');

figure(2);
plot(tick,b,'b-');
hold on;
plot(tick,c,'r+-','MarkerSize',3);
xlabel('\theta_s');
ylabel('RCS(dBsm)');
legend('Original PO-MM','PO-MM with EPA');
title('VV polarization');